Mathematical simulation of the near-bottom section of an ascending twisting flow

The available results of laboratory experiments on the formation of free vortices and controlling of their behavior are compared with the results of mathematical simulation of corresponding flows. This is accomplished by constructing solutions for a set of gas dynamics equations. The comparison is performed for a specific scheme of origination and functioning of free ascending twisting flows. In particular, it is shown that the experimental results confirm the proposed scheme of the origination and initial twisting of ascending vortex flows and validate the reason of their stable functioning with the help of the method intended for controlling generated vortices using vertical grids which was implemented in the experiments. The fact of origination of an ascending flow twisting and its directing is mathematically substantiated using the solution to a specific initially edge problem for a set of gas dynamics equations. A stationary flow whose parameters are close to gas-dynamic parameters of free vortices reproduced in the experiments is calculated. © 2013 Pleiades Publishing, Ltd

Strong magnetohydrodynamic turbulence with cross helicity

Magnetohydrodynamics (MHD) provides the simplest description of magnetic plasma turbulence in a variety of astrophysical and laboratory systems. MHD turbulence with nonzero cross helicity is often called imbalanced, as it implies that the energies of Alfv\'en fluctuations propagating parallel and anti-parallel the background field are not equal. Recent analytical and numerical studies have revealed that at every scale, MHD turbulence consists of regions of positive and negative cross helicity, indicating that such turbulence is inherently locally imbalanced. In this paper, results from high resolution numerical simulations of steady-state incompressible MHD turbulence, with and without cross helicity are presented. It is argued that the inertial range scaling of the energy spectra (E^+ and E^-) of fluctuations moving in opposite directions is independent of the amount of cross-helicity. When cross helicity is nonzero, E^+ and E^- maintain the same scaling, but have differing amplitudes depending on the amount of cross-helicity.Comment: To appear in Physics of Plasma

Collapse Dynamics of a Homopolymer: Theory and Simulation

We present a scaling theory describing the collapse of a homopolymer chain in poor solvent. At time t after the beginning of the collapse, the original Gaussian chain of length N is streamlined to form N/g segments of length R(t), each containing g ~ t monomers. These segments are statistical quantities representing cylinders of length R ~ t^{1/2} and diameter d ~ t^{1/4}, but structured out of stretched arrays of spherical globules. This prescription incorporates the capillary instability. We compare the time-dependent structure factor derived for our theory with that obtained from ultra-large-scale molecular dynamics simulation with explicit solvent. This is the first time such a detailed comparison of theoretical and simulation predictions of collapsing chain structure has been attempted. The favorable agreement between the theoretical and computed structure factors supports the picture of the coarse-graining process during polymer collapse.Comment: 4 pages, 3 figure

Why polymer chains in a melt are not random walks

A cornerstone of modern polymer physics is the `Flory ideality hypothesis' which states that a chain in a polymer melt adopts `ideal' random-walk-like conformations. Here we revisit theoretically and numerically this pivotal assumption and demonstrate that there are noticeable deviations from ideality. The deviations come from the interplay of chain connectivity and the incompressibility of the melt, leading to an effective repulsion between chain segments of all sizes $s$. The amplitude of this repulsion increases with decreasing $s$ where chain segments become more and more swollen. We illustrate this swelling by an analysis of the form factor $F(q)$, i.e. the scattered intensity at wavevector $q$ resulting from intramolecular interferences of a chain. A `Kratky plot' of $q^2F(q)$ {\em vs.} $q$ does not exhibit the plateau for intermediate wavevectors characteristic of ideal chains. One rather finds a conspicuous depression of the plateau, $\delta(F^{-1}(q)) = |q|^3/32\rho$, which increases with $q$ and only depends on the monomer density $\rho$.Comment: 4 pages, 4 figures, EPL, accepted January 200

Long Range Bond-Bond Correlations in Dense Polymer Solutions

The scaling of the bond-bond correlation function $C(s)$ along linear polymer chains is investigated with respect to the curvilinear distance, $s$, along the flexible chain and the monomer density, $\rho$, via Monte Carlo and molecular dynamics simulations. % Surprisingly, the correlations in dense three dimensional solutions are found to decay with a power law $C(s) \sim s^{-\omega}$ with $\omega=3/2$ and the exponential behavior commonly assumed is clearly ruled out for long chains. % In semidilute solutions, the density dependent scaling of $C(s) \approx g^{-\omega_0} (s/g)^{-\omega}$ with $\omega_0=2-2\nu=0.824$ ($\nu=0.588$ being Flory's exponent) is set by the number of monomers $g(\rho)$ contained in an excluded volume blob of size $\xi$. % Our computational findings compare well with simple scaling arguments and perturbation calculation. The power-law behavior is due to self-interactions of chains on distances $s \gg g$ caused by the connectivity of chains and the incompressibility of the melt. %Comment: 4 pages, 4 figure

Single Chain Force Spectroscopy: Sequence Dependence

We study the elastic properties of a single A/B copolymer chain with a specific sequence. We predict a rich structure in the force extension relations which can be addressed to the sequence. The variational method is introduced to probe local minima on the path of stretching and releasing. At given force, we find multiple configurations which are separated by energy barriers. A collapsed globular configuration consists of several domains which unravel cooperatively. Upon stretching, unfolding path shows stepwise pattern corresponding to the unfolding of each domain. While releasing, several cores can be created simultaneously in the middle of the chain resulting in a different path of collapse.Comment: 6 pages 3 figure

Logarithmic corrections of the avalanche distributions of sandpile models at the upper critical dimension

We study numerically the dynamical properties of the BTW model on a square lattice for various dimensions. The aim of this investigation is to determine the value of the upper critical dimension where the avalanche distributions are characterized by the mean-field exponents. Our results are consistent with the assumption that the scaling behavior of the four-dimensional BTW model is characterized by the mean-field exponents with additional logarithmic corrections. We benefit in our analysis from the exact solution of the directed BTW model at the upper critical dimension which allows to derive how logarithmic corrections affect the scaling behavior at the upper critical dimension. Similar logarithmic corrections forms fit the numerical data for the four-dimensional BTW model, strongly suggesting that the value of the upper critical dimension is four.Comment: 8 pages, including 9 figures, accepted for publication in Phys. Rev.

The Bak-Tang-Wiesenfeld sandpile model around the upper critical dimension

We consider the Bak-Tang-Wiesenfeld sandpile model on square lattices in different dimensions (D>=6). A finite size scaling analysis of the avalanche probability distributions yields the values of the distribution exponents, the dynamical exponent, and the dimension of the avalanches. Above the upper critical dimension D_u=4 the exponents equal the known mean field values. An analysis of the area probability distributions indicates that the avalanches are fractal above the critical dimension.Comment: 7 pages, including 9 figures, accepted for publication in Physical Review

Vortex Entanglement and Broken Symmetry

Based on the London approximation, we investigate numerically the stability of the elementary configurations of entanglement, the twisted-pair and the twisted-triplet, in the vortex-lattice and -liquid phases. We find that, except for the dilute limit, the twisted-pair is unstable and hence irrelevant in the discussion of entanglement. In the lattice phase the twisted-triplet constitutes a metastable, confined configuration of high energy. Loss of lattice symmetry upon melting leads to deconfinement and the twisted-triplet turns into a low-energy helical configuration.Comment: 4 pages, RevTex, 2 figures on reques

Stochastic Flux-Freezing and Magnetic Dynamo

We argue that magnetic flux-conservation in turbulent plasmas at high magnetic Reynolds numbers neither holds in the conventional sense nor is entirely broken, but instead is valid in a novel statistical sense associated to the "spontaneous stochasticity" of Lagrangian particle tra jectories. The latter phenomenon is due to the explosive separation of particles undergoing turbulent Richardson diffusion, which leads to a breakdown of Laplacian determinism for classical dynamics. We discuss empirical evidence for spontaneous stochasticity, including our own new numerical results. We then use a Lagrangian path-integral approach to establish stochastic flux-freezing for resistive hydromagnetic equations and to argue, based on the properties of Richardson diffusion, that flux-conservation must remain stochastic at infinite magnetic Reynolds number. As an important application of these results we consider the kinematic, fluctuation dynamo in non-helical, incompressible turbulence at unit magnetic Prandtl number. We present results on the Lagrangian dynamo mechanisms by a stochastic particle method which demonstrate a strong similarity between the Pr = 1 and Pr = 0 dynamos. Stochasticity of field-line motion is an essential ingredient of both. We finally consider briefly some consequences for nonlinear MHD turbulence, dynamo and reconnectionComment: 29 pages, 10 figure